Exact and Asymptotic Conditions on Traveling Wave Solutions of the Navier-Stokes 

Equations 



Y. Charles Li 

Department of Mathematics, University of Missouri, Columbia, MO 65211, f/S^Q 



o 
o 

(N 



Divakar Viswanath 

Department of Mathematics, University of Michigan, Ann Arbor, MI 48109, US$\ 

(Dated:) 

We derive necessary conditions that traveling wave solutions of the Navier-Stokes equations must 
satisfy in the pipe, Couette, and channel flow geometries. Some conditions are exact and must hold 
for any traveling wave solution irrespective of the Reynolds number (Re). Other conditions are 
asymptotic in the limit Re — » oo. The exact conditions are likely to be useful tools in the study of 
transitional structures. For the pipe flow geometry, we give computations up to Re = 100000 showing 
the connection of our asymptotic conditions to critical layers that accompany vortex structures at 
high Re. 
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Turbulence is sometimes said to be the last unsolved 
problem of classical physics. However, in a sense it is a 
fully solved problem, since we know with near certainty 
that the Navier-Stokes equations (NSE), along with the 
no-slip boundary conditions [l[ , are an excellent physical 
model for all the phenomena associated with turbulence 
and transition. 

Although the physical model is known, figuring out 
the solutions of the NSE has proved to be a very difficult 
problem, computationally or analytically. To illustrate 
this difficulty, we note that while useful models of the 
structural strength of an entire aircraft can be built using 
1 to 100 million elements, it takes more than a billion 
grid points to completely resolve the flow in a cubical 
volume of extent an inch 3 at the surface of a car moving 
at around 60 m.p.h. 

The Navier-Stokes equations, which model the evolu- 
tion of the velocity field u of an incompressible fluid, are 
du/dt + (u.V)u = — Vp + (l/i?e)Au, where the veloc- 
ity field u must satisfy the incompressibility constraint 
V.u = 0. There is no explicit equation for evolving the 
pressure p, and Re denotes the Reynolds number. 

Traveling wave solutions of the form u(x, t) = u(x — ct) 
form our main topic. As the motion of turbulent fluids 
is characterized by disordered and intermittent fluctua- 
tions about a mean, the significance of traveling wave 
solutions may seem limited. Indeed, one has to look at 
more complicated solutions to begin to understand tur- 
bulent fluctuations [j|. However, there is gathering evi- 
dence that traveling wave solutions may help understand 
certain coherent structures in transitional pipe flow 

Further, using certain lower-branch traveling waves, we 
can exhibit critical layers at high Re [H, 0] that are far 
beyond the reach of ordinary direct numerical simula- 
tion. The Orr-Sommerfeld equation, which governs the 
propagation of infinitesimal normal mode perturbations 
of a base flow, is singular at Re = oo, which is the origin 
of the theory of critical layers [l(J 11 1. The ability to 
compute critical layers in fully resolved numerical solu- 



tions of the NSE could be significant, as critical layers 
occur in many important situations [llfl . The gigantic 
trailing vortices that escape from the boundary layers of 
airplanes during take-off may develop critical layers [l2| . 
So could vortices shed by wind turbines, with possible 
implications for the optimal arrangement of turbines in 
a wind farm. 

Many linearly unstable (and non-laminar) traveling 
wave solutions and equilibria (which are special cases 
with c = 0) of Couette [H [jj], channel [ll, HI], and 
pipe [13, EH geometries have been computed. A notably 
systematic and extensive effort is due to Gibson and oth- 
ers [Til HH- The exact conditions we derive must hold 
for the velocity fields of all these solutions. In addition, 
analogous conditions must hold for periodic solutions and 
relative periodic solutions that travel only in the stream- 
wise direction. 

The mere existence of traveling wave solutions does 
not imply their relevance to phenomena as they manifest 
themselves in Nature and in technology. However, efforts 
to connect traveling waves computed in short pipes to 
puffs have been partially successful 0, 0] ■ Puffs are tran- 
sitional structures observed in pipes around Re = 2000 
that are approximately 20 pipe diameters in axial length 
and which travel with a speed that is somewhat less than 
the mean streamwise velocity [2l|. Tantalizingly, there 
are hints that an entire puff may correspond to a trav- 
eling wave or some such solution of pipe flow [22l | . Our 
exact conditions will be helpful in the investigation of 
that possibility. 

We now turn to the derivation of the exact conditions. 
In the velocity field u = (u,v,w) of the NSE, u, v and 
w are the streamwise (coordinate axis x), wall- normal 
(y), and spanwisc (z) components, respectively, for the 
rectangular Couette and channel geometries. In the case 
of pipe flow, u, v and w are the radial (r), polar (6) and 
streamwise (or axial) (z) components, respectively. 

In plane Couette flow, the walls at y = ±1 move in 
the x direction with speeds equal to ±1. The boundary 
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conditions in the streamwise and spanwise directions are 
periodic, with the periods taken to be 2ith x and 2ith Zl re- 
spectively. For pipe flow, we assume the axial or stream- 
wise boundary condition to be periodic with period 2itA. 
The walls are no slip in all cases. The derivations are 
given mainly for the plane Couette flow geometry. 

Let u(x, t) = u(x — ct) be a traveling wave solution of 
plane Couette flow. We assume c = (c, 0,0) so that the 
traveling wave moves in the streamwise direction only. If 
u = (u, v, w), the streamwise or x component of the NSE 
gives 



cd x u + (ud x u + vdyU + wd z u) 



in 

Here c p /Re gives the pressure gradient in the stream- 
wise directions, with c p — for plane Couette flow 
and c p > for channel flow. Let U (y, z) denote 
(2ttA x )~ 1 Jq^ A * u(x,y, z) dx, which is the mean stream- 
wise component of u. From |T]), we get 



(ud x u + vdyU + wd z u) 



AU 



Re 



(2) 



where the overline denotes streamwise averaging. At the 
walls y = ±1, d x u = and v = w = because of no-slip. 
For the same reason, d zz U = at the walls. Therefore, 



Cp + OyyU = 



(3) 



must hold at the walls. 

As the velocity field u has zero divergence, we may 
rewrite ([2]) as 



V • (u 2 , uv, uw) = dyuv + d z uw 



AU 



Re 



(4) 



If (j4j) is integrated over the cross-section, Green's theo- 
rem applies to the expression in the middle of jl]). The 
integral of the middle term must be zero because v = 
at the walls and uw is periodic in z. Thus we have 



2ttA 2 r + l 



( AU + c p ) dy dz = 0. 



(5) 



The derivation of the necessary conditions ([3| and (|5|) 
applies to channel flow with no change. However, c p ^ 
for channel flow. 

The conditions Q and §5§ must be satisfied by all trav- 
eling wave solutions of plane Couette flow or channel flow, 
whose wave speed vector c only has a streamwise compo- 
nent. Indeed, those conditions must be satisfied by all pe- 
riodic solutions u(x, t) = u(x, t + T), T being the period, 
or relative periodic solutions u(x, t) = u(x + s, t + T) if 
the shift s only has a streamwise component. To form 
U in those instances, one must average both over a sin- 
gle period and in the streamwise direction as a simple 
modification of our derivation will show. 



For the case of pipe flow, let c = (0, 0, c) so that the 
traveling wave travels in the streamwise direction only. 
Let W(r,B) = (27rA)- 1 J^ 1 w(r,9,z)dz be the mean 
streamwise velocity. The analogue of ([3]) requires 



d r W 



= 



(G) 



at all points on the circumference. If we assume the pipe 
radius to be 1, the analogue of (O is 



(AW ■ 



Cp)r dr d9 



0. 



(T) 



The derivation of the necessary conditions ((6|) and j7|) for 
pipe flow is similar to that of their Couette analogues. 

Traveling waves normally arise from saddle-node bi- 
furcations with increasing Re [l3|, 15, 17, 18, H^]. The 



branch corresponding to lower energy dissipation is called 
the lower branch. We will now derive certain scalings 
with respect to increasing Re that are characteristic of 
the lower branch families. 

In the case of plane Couette flow or channel flow, if a 
traveling wave solution is given by u(x — ct), the velocity 
field u(x) can be decomposed as 



uo(y,z) + ^2(u n (y,z)exp(iax) 



(8) 



where a — l/A x . We take Uo = (U,vq,wo) and 
Uj = (ui,Vi,Wi) for i > 1. For pipe flow, the de- 
composition analogous to (J8]) is given by Uo(r, 6) + 
J2^Li{ u n(r,0) exp(iaz) +c.c), with a = 1/A. We take 
u, = (ui,Vi,Wi) for i > 1 as for Couette flow, but 
Uo = (uo, vq, W) for pipe flow. 

The scalings of the lower branch family that are known 
or that will be derived apply to the mean streamwise 
velocity (U or W), or the rolls ((vo,wq) or (u ,v j), or 
the magnitude of modes such as ui. It is an empirical 
fact (but see fl8L |25| ) that the rolls and the Ui mode 
diminish in magnitude approximately at the rate i?e _1 . 
Higher modes with n > 1 appear to diminish even faster. 
The derivations assume these scalings. In addition, the 
dissipation rate of the lower branch families decrease with 
increasing Re, assuming that the dissipation rate of the 
laminar solution is normalized to be 1 [1, Q . 

The wall-normal or y component of the n = 1 mode of 
the NSE gives 

ia(U-c)vx = -d y p! + Re~ 1 (-a 2 v 1 +dyV 1 +d z v 1 )-\ 

(9) 

for plane Couette or channel flow. The first neglected 
terms in ([9]) are —v$d y vi — wq8 z vi — v\d y Vo — wid z vo- 
The analogous equation for the radial component of pipe 
flow is ia(W — c)u\ = —d r pi + Re^ 1 A r u\, where A r 
corresponds to the usual form of the Laplacian in the 
radial component of the NSE. Terms such as —vxdgvo/r 
are neglected. 
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Re = 2600 



Re = 1 0000 



Re = 1 00000 




FIG. 1: Three traveling wave solutions of pipe flow are shown. Units are chosen to make the pipe radius and the centerline 
speed of the laminar solution equal to 1. The boundary condition fixes the mass flux (and not the pressure gradient) to be 
that of the laminar solution. Solutions (a) and (b) correspond to the c family and the upper dash-dot curve of Figure 7 of [l8| ]. 
respectively. For more on (c), see In the upper plots, the four contour levels of |ui| are equispaced in (0,max), with max 
being .0094, .0011, and .0002, respectively. The thick line is the critical curve. The plots below correspond to the ones directly 
above. Plot (d) (the thick line is for c p ) can be validated against Figure 7 of ijij after replacing Re and c by c p Re and c/c p , 
respectively. In (e) and (f), the legends show the choice of f(x) and the curves illustrate (| 14f) (a ~ 1.5). 



Using (|9|), Walcffc ct al.Q estimated that most of the 
variation in v\ is concentrated in a region around the 
critical curve U = c, with the width of that region scaling 
as i?e -1 / 3 . In the case of pipe flow, an identical argument 
gives W = c as the equation of the critical curve. The 
top set of plots in Figure [1] illustrate the critical layer in 
the case of pipe flow. 

To derive further scalings, we consider the streamwise 
component of the n = mode of NSE, which is 



vodyU + wod z U 



Rer 1 (c p + AU) +M (10) 



for plane Couctte or channel flow. The pipe flow analogue 
is u d r W + (v /r)d e W = Re~ x {cp + AW) + M. In dTU]) 
and its pipe flow analogue, — M equals the n = mode 
of the streamwise component of ((u — Uo).V)(u — Uo). 
From here on we restrict the derivation to plane Couette 
flow or to channel flow. 

Since Uq = (U, vq, wq) has zero divergence, we can find 
a function ?/>(?/, z) such that Vq = d z ip and wq = —dyij). 
We then get 



where 



Ltp = Re-\c p + AC/) + 71/, 



(d v u) d z - (d z u) d y 



(ii) 



(12) 



The skew-symmctry of the linear operator L is the key 
to deducing further scalings. The skew-symmetry of L is 
likely to be important in attempts to find an asymptotic 
theory for the critical layer. 

Let <p(y, z) and ij)(y, z) have z periods of 2t:A z and 
be sufficiently smooth. The following calculation uses 
integration by parts: 



i2ttA 



4>i/jU z \ dz 



2ttA pi rl 

<j)Ltj) dy dz = I (pipUy^'^dy 
l J-i 

/ (ipU y <f) z + i/j(j)Uy Z )dy dz - 
o J-i Jo 

2ttA p\ 

j {i>U z (j>y + ij)(j)Uy Z )dydz, 
o J-i 



where the subscripts are for partial derivatives. On the 
right hand side, two double integral terms cancel and the 
single integral terms are both zero because U y is periodic 
in z and U z is zero at the walls (from no-slip). We are 
left with — J J ipLtfidydz on the right, verifying skew- 
symmetry of the operator L. 

From direct substitution into (fl2|) . it is evident that 
L(f(U)) = for any smooth /. Thus the functions f(U) 
are in the kernel of the anti-symmetric operator L. Since 
the linear system (|11[) can be solved for ip (or equivalcntly 
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for the rolls), the Fredholm alternative implies that 
J J f(U)(cp + AU + ReM)dydz = 0. 

For lower-branch traveling wave families with u — Uo of 
magnitude Re~ a with awl, the magnitude of M is 
approximately Re~ 2 in the limit Re — ► oo. We have 

J f(U)(AU{y,z) + c^dy dz = 0(Re~ a ) (13) 

for any smooth /. Here a > 1 is possible if there are 
cancellations in the integral of M over the cross-section. 
The analogous condition for pipe flow is given by 

J f(W) (AW(r, 9) + Cp) rdr d6 = 0(Re~ a ) (14) 

with a as above. 

In addition to the pipe families of Figure [T] we com- 
puted a lower branch equilibrium family and a traveling 
wave family up to Re = 45000 and Re = 7000. The 
(a) and (b) families of Figure [1] could not be continued 
to Re much higher than shown in the top plots. For a 
given resolution, we cannot expect to find solutions if 
the rolls, which diminish in magnitude with Re, are too 
small to be detected. Even after using sufficient resolu- 
tion, the GMRES-hookstep iterations (see @,@|) became 
very slow. Even though the residual error could be made 
quite small, the norm of the Newton steps became quite 
large and increased with iteration. Although it is uncer- 
tain if the lower branch families exist in the Re —> oo 
limit, Figure [1] amply demonstrates that they exist for 
large enough Re for the predicted scalings to hold. 

The critical curves are away from the pipe walls and 
have an inward indentation where the counter-rotating 
vortices face each other. Since the critical behavior is ev- 
ident even for Re = 2600, we suspect that critical curves 
or surfaces may give a way to visualize the structure of 
puffs in transitional pipe flow. 

In summary, we have given a number of necessary con- 
ditions for equilibrium, traveling wave, periodic, and rel- 
ative periodic solutions in plane Couctte, channel and 
pipe geometries. We have argued for the importance of 
critical layers in high Re fluid flow and shown the con- 
nection of our analysis to critical layers. In addition, the 
conditions that we have derived are likely to be useful in 
the study of transitional structures such as puffs in pipe 
flow. 
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